Spin nematic states in antiferromagnets containing ferromagnetic bonds 
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The majority of recent works devoted to spin nematic phases deal with either frustrated magnets 
or with those described by Hamiltonians with large non-Heisenberg terms. We show in the present 
study that non-frustrated antiferromagnets (AFs) containing ferromagnetic (FM) bonds can show 
nematic phases in strong magnetic field. Among particular spin systems discussed are a ladder with 
FM rungs, two AF layers coupled ferromagnetically, a chain containing alternating AF and FM 
bonds and an AF anisotropic spin-1 chain. 
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Introduction. — Frustrated spin systems have offered in a) 

recent years a wealth of opportunities for the study of a 
broad range of novel types of states and phase transitions. 
Spin nematic phases form a class of objects in this area 
which has received much attention. Spin nematic states 
are spin-liquid-like states which show a multiple-spin or- 
dering without the conventional long-range magnetic or- 
der. The two-spin ordering can be generally described 
by the tenso^ Q?f = (S?S?) - S aj3 (SjSi)/3. The anti- 
symmetric part of Q^f is related to the vector chirality 
(Sj x Si) and describes a vector chiral spin liquid which c \ 

can be stabilized in quantum spin models at T = by 
a sizable ring-exchange, 2 and can be found in classical 
frustrated spin systems at T ^ OP The symmetric part 



of Q**i describes a quadrupolar order which has been ex- 
tensively studied both theoretically and experimentally 
in frustrated systems with ferromagnetic (FM) and an- 
tiferromagnetic (AF) nearest-neighbor and next-nearest- 
neighbor couplings, respectively, in strong magnetic field 
H (see, e.g., Ref. 4 and references therein) and in mag- 
nets with large non-Heisenberg spin couplings such as 
biquadratic exchange (S1S2) 2 . 5 It has been also shown 
recently that quantum fluctuations accompanied by a siz- 
able single-ion easy-axis anisotropy can also stabilize a 
nematic phase in the kagome spin-1 antiferromagnet. 6 
This study was motivated by recent experiments on 
Ni 3 V 2 8 . 

It is well established that the attraction between 
magnons caused by frustration is the origin of quadrupo- 
lar and multipolar phases in quantum magnets P In par- 
ticular, the bottom of the one-magnon band lies above 
the lowest multi-magnon bound state at H = H S1 where 
H s is the saturation field, as a result of this attraction 
in magnets with FM and AF couplings between nearest- 
and next-nearest neighbors, respectively. Then, transi- 
tions to nematic phases at H < H s in such systems are 
characterized by a softening of the multi-magnon bound- 
state spectrum rather than the one-magnon spectrum. 

We show in the present paper that FM bonds in 
AF non-frustrated spin systems can also lead to the 
magnon attraction and to stabilization of nematic phases 
in strong magnetic field. To be specific, we discuss spin- 
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FIG. 1. (Color online.) (a)-(c) Spin-| systems which show 
a nematic phase in strong magnetic field and which are dis- 
cussed in the present paper in some detail. Antiferromagnetic 
and ferromagnetic bonds are shown in black and blue, respec- 
tively, (d) A sketch of the phase diagram at small T of the 
considered systems. The canted or the nematic phase is ab- 
sent in some cases (see the text). A small inter-chain, inter- 
bilayer and inter- ladder couplings are taken into account. 



I ladder with FM rungs, two AF layers coupled ferro- 
magnetically and a chain containing alternating AF and 
FM bonds which are presented in Fig. [lja)-(c). The 
spin ladder with FM rungs and the alternating chain 
have received considerable interest recently which has not 
been related, however, to their nematic behavior in mag- 
netic field (see Refs. ^ and RefsP, respectively, and ref- 
erences therein). Although interaction between spins is 
antiferromagnetic in the majority of compounds contain- 
ing weakly coupled spin chains, ladders or bilayers, some 
(mainly organic-based) materials containing FM bonds 
have been synthesized recently^ that has stimulated the 
theoretical activity in this field. 

Model and technique. — All spin-^ systems under dis- 
cussion are described by the Hamiltonian 

T-i = J^2 (Si.zSi.z + S 2 ,iS 2 ,0 - j ^2 s M s 2,i 

(hi) i 

- a Y, sii s ii - H H ( S U + S W + «'> (i) 



where S n ^ is the n-th spin in the i-th FM bond, j and 
A are values of the isotropic and anisotropic parts of the 
FM exchange coupling, respectively, (i, I) denote nearest 
neighbor spins coupled antiferromagnetically, we set the 
AF coupling constant J to be equal to unity in our calcu- 
lations, and %' describes small inter-ladder (inter-bilayer 
or inter-chain) interaction. 

The ground state of Hamiltonian (fil) has a collinear 
antiferromagnetic spin structure at small H and the 
magnon spectrum has a gap induced by the easy-axis 
anisotropy A. There is a sequence of phase transitions 
upon the field increasing. 

If the easy-axis anisotropy A is smaller than a crit- 
ical value A c (the classical value of A c is Jo, where 
J p = J V • exp(zpR j7 ), that does not depend on j), the 
first phase transition is the first-order spin-flop one that 
happens when the field reaches a value H s f. One has for 
the classical value of H s f 



H? f = ±y/A(2J -A). 



(2) 



There is a canted AF spin structure at H > H s f and 
one of the magnon branches is gapless (as a consequence 
of the continuous SO (2) symmetry breakdown in this 
phase). Our experience suggests that there would be 
only one phase transition to the fully polarized phase 
at H = H s upon further field increasing. However, we 
show below that in a range of parameters the fully po- 
larized phase is preceded by a nematic one with the or- 
der parameter (S^S^) (see Fig. lfd)). Consideration 
of the nematic order parameter symmetry shows that 
SO(2)/Z 2 symmetry is broken in the nematic phase (see, 
e.g., RefP^). Then, there is a Goldstone mode in the ne- 
matic phase and the transition from the canted AF phase 
to the nematic one is apparently of the 2D Ising type be- 
cause the Z2 subgroup breaks down. 

If A is large enough (e.g., at A ^> J, j), there is no 
canted phase. We show below that even in this case the 
collinear and the fully polarized gapped phases are sepa- 
rated by the gapless nematic one. 

We examine in the present paper the possibility of the 
nematic phase formation just below the saturation field 
H s by considering the transition from the fully polarized 
phase. We use for this purpose the approach suggested 
recently by one of us in Ref.^for high- field nematic phases 
analysis that is based on the following bond-operator spin 
representation: 



S l,i =b\(Li +Ci, 
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Kb, 



where aj, b\ and c\ are Bose-operators which create 
three spin states from the vacuum |0) = |tt) as follows: 

<4l°) = \U), b\\0) = |t;), and cj|0) = ||t), where all 
spins have the maximum projection on the field direc- 
tion at the state |0). To avoid contribution of unphysical 



states containing more than one particle a, b or c on a 
FM bond, we add to the Hamiltonian constraint terms 
describing an infinite repulsion between particles on each 



FMbond: U E*K a J a ^HA T W 



uthh 



-c]c]ciCi+a]b]aibi 



b]c]biCi + a]c]aiCi), where U — >■ +00. In particular, it is 
shown in Ref J^ that (i) this approach is quite convenient 
for discussion of the quantum phase transition to the ne- 
matic phase, (ii) along with some new results it yields 
also those obtained by other methods P^U an d (iii) al- 
though it gives quantitatively correct results at H « H s 
when H < H s , it works in the nematic phase qualita- 
tively also when H is not very close to H s . For the sake 
of self-consistency, we present below some detail of this 
approach. We confirm below our key analytical results 
by numerical ones obtained using finite cluster diagonal- 
ization technique.^ 

Substituting Eqs. pi into Hamiltonian (fTJ, taking into 
account the constraint terms and neglecting for a moment 
%' one obtains for the ladder and the bilayer 

^ = E (( 6 p 6 p + c p c p) ( H + ^( J p - J °) + \u + A ) 



- 2"i (&pCp + 4> b v) + «p a p( 2 # " J o)J , (4) 

U 3 = —7= V (a\b 2 c 3 (J 2 + J 3 ) + b\4a 3 (Jj + J 2 )) , (5) 

2v N — V ' 

% 4 = — ^2 ( a 14 a 3«4 ( J1-3 + U) + Ub\clb 3 cs 

+ \a\bla3b4 + a\cla 3 cA I J1-3 + n Jl ~ 4 + U ) 

+ (b\bib s b 4 + 44C3C4) Q J1-3 + u)\ (6) 



where N is the number of FM bonds, the momentum 
conservation laws J2i Pi = are implied in Eqs. J5| and 
d6J), and we omit some indexes p in Eqs. ^ and mf. It is 
convenient to introduce the following Green's functions: 



G a {k) = 


-i{a k a\), 






(7) 


G b (k) = 


-*(&*&£>> 


G c (k) = 


-i(ckc{), 


(8) 


F(k) = 


-i(b k c{), 


F(k) = ■ 


-i(c k b{). 


(9) 



Poles of G a {k) give the spectrum of a particles, whereas 
poles of Gb(k), G c (k), F(k) and F(k) which have the 
same denominator determine spectra of two one-magnon 
branches.^ 

Particles b and c carrying spin 1 are of the one-magnon 
nature. Their spectra calculated at H > H s using 
Eqs. (|3| coincide with those derived using the conven- 
tional approaches such as the Holstein-Primakoff trans- 
formation. It can be shown 4 that spectra of one-magnon 
excitations are determined solely by H2 (Eq. Q) at 
H > H s and they have the form 



e 1 (p) = J ff+-(J p -J ) 
£2(p) = ei(p)+j. 



:A, 



(10) 
(11) 
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FIG. 2. (Color online.) a) Diagrams for the self-energy part 
S (A;) of a-particles at H > if s . Green's functions of b- and c- 
particles Gb(p), G c (p), F(p) and F(p) are denned by Eqs. (J8|- 
([9}. Circles represent renormalized vertices. Equations for the 
vertices are presented in panel (b). Bare vertices are defined 
by Eqs. (5| and Q. 



The lower branch ei(p) has a minimum at p = tt (or 
(7r, 7r)) and a gap which closes when H becomes equal to 



H C = J - A/2. 



(12) 



In contrast to b and c particles, a particles carrying 
spin 2 are of the two-magnon nature. Their spectrum 
coincides with the two-magnon bound-state spectrum 
found using conventional methods!^ To find the spec- 
trum e a (p) of a particles at H > H s one has to take 
into account diagrams shown in Fig. [2|a) which contain 
three-particle vertexes. Equations for them are presented 
in Fig. f2^b). The minimum of e a (p) is at p = and 
the gap in e a (p) closes at H = H' c . If H c > H r c , the 
transition takes place at H = H c = H s to the canted 
phase which can be described as the Bose-Einstein con- 
densation (BEC) of one-magnon excitations. As it is seen 
from Eqs. pi), (S^) becomes finite in this case, where 
_L denotes the projection on the xy plane. In contrast, 
if H c < H' c the transition from the fully polarized phase 
to the nematic one takes place at H = H' c = H s which 
can be described within our approach as the BEC of a 
particles. It is seen from Eqs. (J3| that (S^) = in the 

nematic phase and (S^jS^j) = (a]) ex e~ i(t> ^H s (T) -H 
is the nematic order parameter (for H! 7^ 0), where <j) is 
an arbitrary phase. 4 All the static two-spin correlators 
decay exponentially in the nematic phase. We calculate 
below e a (p) for the selected spin systems and find stabil- 
ity conditions of the nematic states. 

AF ladder with FM rungs (see Fig. \j\a)). — 
Unfortunately, the general expression for e a (p) is quite 
complicated but it is simplified greatly in the limiting 
case of the Ising exchange on FM bond (i.e., at j — 0). 
To illustrate main properties of the nematic phase we 
consider here in some detail the case of j = 0. In partic- 
ular, one finds at %' — and H > H s 



e a (p) = 2H-2- 



4 cos 2 



■A 2 . 



(13) 



The gap in e (p) vanishes at H s = 1 + y/l + A 2 ~Jl- A/2 
that is larger than H c = 2 — A/2 (see Eq. (12)) at all 



positive A. Then, the nematic phase is stable at A > 
when j = 0. Interestingly, the nematic phase arises also 
in the case of A ^> 1, when the canted phase is absent. 
Thus, the gapless nematic phase separates two gapped 
phases at large A, the fully polarized and the collinear 
ones. 



Spectrum ([13]) has a quadratic dispersion near its min 

V4 



imum at p = 0: e a (p) 
F>\\p 2 , where D 



2H - 2 + A - V4 + A 2 + 
A 2 . The spectrum remains 



1/2V4 
quadratic also after taking into account a small %' that 
leads to the 3D BEC relation H s (0) - H S (T) ex T 3 / 2 P1 

One obtains in the first order in p and in the leading 
order in the inter-ladder interaction as it is done in Ref. 4 
for quasi- ID frustrated systems 



\^l,i^2,j+n) 



\ S l,j S lj+n) - \ S 2,j S 2j+n) 




(14) 




where n > and p = {a\ai) ex (H S (T) — H) is the "con- 
densate" density. The rest two-spin static correlators 
containing S + or £" are equal to zero and (S'f^-S'f - +n ) = 



2,j^2,j+n) 



p . It is seen from Eqs. (14)-(16) that all 



correlators decay exponentially at A > 0. 
One obtains for the magnetization 



Wj) 



8 16\/4 + A 2 



A 2 



A 5 



(17) 



where q = 1,2. It may seem that Eqs. (14)-(17) are 



invalid for arbitrary small A. However, the above results 

are valid in the near vicinity of H s (i.e., at very small p) 

in this case because H s and H c merge in the limit A — >• 0. 

It is implied in Eqs. (Il4|-([l7|) that W ^ so that 



p 7^ 0. The situation is completely different in the purely 
ID case because p = and there is no long range ne- 
matic order. Bearing in mind the quadratic dispersion of 
e a (p) near its minimum and using results of ID Bose-gas 
discussions 14 one obtains 13 at H « H* and T = 



H^H^-K 1 /2(//s U) 



(S z q , J+n (t)S^(0)) « (5?) 
1/1 1 




(18) 
(19) 



7r\(n + iut) 2 (n — iut) 2 

(^oWs+ (t)s M (o)s 2 ; n (o)) 



B 



cos (n (a jCLj)n) 
n 2 + u 2 t 2 



Bo 



where q, p = 1, 2, n — >• 00, u = 47rD|| (aja^) and B12 are 
constants. 



y/\n + iut\ 
i)aj) 



(20) 




FIG. 3. (Color online.) Regions are shown of the nematic 
phase stability in the ladder with FM rungs (inequality (21 )), 
in the bilayer with FM coupling (inequality (22)), and in the 
alternating chain. 



Analysis of the general expression for e a (p) gives the 
following general criterion of the nematic phase stability: 



A> 



4j 



j + 2v7(JT2) 



(21) 



that is shown graphically in Fig. [3J One concludes from 
Eq. (21) that only large A > j can stabilize the nematic 
1. In contrast, quite small anisotropy on FM 
1 is sufficient for this purpose if j ^> 1. 



phase if j 
bonds j ^> A 

In the limiting case of j — »> oo, our model describes the 
spin-1 AF chain with the easy- axis single-ion anisotropy 
A. The nematic order parameter reads in this case as 
((S^) 2 ) and one obtains from Eq. (21) that the nematic 
phase is stable at A > 4/3. Similar nematic phase is 
discussed in RefPin spin-1 kagome AF with large single- 
ion easy- axis anisotropy at H = 0. 

Our finite cluster calculations confirm that the transi- 
tion takes place at H — H s to the nematic phase when 
inequality (21) holds. This numerical consideration is 
simplified by the fact that the Hamiltonian (HI) commutes 
with the z component of the total spin S z and with the 
Zeeman term. As a result all the Hamiltonian eigen- 
states can be classified by eigenvalues M of S z . Let us 



denote E(M) the minimum energy in each M sector at 
H = 0. The ground state energy of a cluster with L 
spins in magnetic field is given by the minimum value of 
E(M) — HM. An important observation is that values 
E(M) — HM at even M sat — M are smaller than those 
with odd M sat — M when H is close to its saturation 
value, where M sat = L/2. Thus, one can expect that a 
condensation takes place in the thermodynamic limit of 
elementary excitations carrying spin 2. Then, the lowest 
state in each even-(M sat — M) sector has zero momentum 
(if the periodic boundary conditions are applied) that is 
also in agreement with the bound states condensation 
scenario. Values of H s obtained numerically as a result of 
analysis of clusters with L = 8-1-48 are in excellent agree- 
ment with the corresponding values found analytically. 
Numerical consideration of clusters with L = 8 -=- 26 simi- 
lar to that performed in Refs.^ confirms also the validity 
of Eq. ( 18 ). For example, one obtains for j = and A = 2 
that \ - (Slj) = a(H s - H) 1 ^ 5 , where a = 1.04 ± 0.11 
and S = 2.05 ± 0.05, that is in good agreement with the 
analytical result (alaj) ~ 1.07 y/H 3 — H. 

AF bilayer with FM coupling (see Fig. \^fb)). — The an- 
alytical analysis of this system is much more complicated 
than that carried out above for the ladder. Then, we re- 
strict ourselves by discussion of the region of the nematic 
phase stability that is defined by inequality (see Fig. |3| 



A> 



m 



Sjf(j) 



(8 + j)/(j)-l' 

■2tt 



1 



dxdy 



(2tt) 2 J j + 2 + cosx + cosy 



(22) 
(23) 



One concludes from Eq. (22) that, similar to the ladder 
with FM rungs, only large A > j can stabilize the ne- 
matic phase if j ~ 1. In contrast, quite small anisotropy 
on FM bonds j ^> A ~ 1 is sufficient for this purpose if 

Alternating chain, (see Fig. wYc)) — Particular expres- 
sions are quite cumbersome in this case. Then, we re- 
strict ourselves by graphical representation of the ne- 
matic phase stability region which is also shown in Fig. [3] 
It is seen that the graphic resembles those for the ladder 
and the bilayer. Numerical consideration confirms exis- 
tence of the nematic phase and the validity of Eq. (18). 
Numerical values of H s are in excellent agreement with 
analytical results. 

To conclude, we demonstrate that FM bonds in non- 
frustrated antiferromagnets can lead to nematic spin 
states in strong magnetic field. The uniaxial anisotropy 
on FM bond is necessary for the nematic phase stabi- 
lization in all the considered systems. The present study 
should stimulate further theoretical and experimental ac- 
tivity in nematic phases discussion both in the considered 
systems and in other ones containing FM bonds. 
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